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Abstract 

We summarize several decades of work in finding values for the percolation threshold pc for site percolation on 
the square lattice, the universal correction-to-scaling exponent Q., and the susceptibility amplitude ratio jC^, in 
two dimensions. Recent studies have yielded the precise values pc - 0.59274602(4), Q. - 72/91 ~ 0.791, and 
C^/C^ = 161.5(2.0), resolving long-standing controversies about the last two quantities and verifying the widely 
used value pc - 0.592746 for the first. 



1. Introduction 

Percolation concerns the formation of long-range connectivity in a system. When the site or bond occupancy 
reaches a critical threshold p^, infinite connectivity first forms. For some systems (e.g., bond percolation on the square 
lattice, or site percolation on the triangular lattice — all in two dimensions (2d)), p^ is known exactly, while for others, 
it must be determined numerically. A notable example of the latter is site percolation on the square lattice, whose 
threshold has been estimated in dozens of studies, many quite extensive, since 1960. In this paper we list some of 
the methods that have been used in making those estimates, and summarize the values that have been found. We also 
summarize results for two related quantities — the correction-to-scaling exponents, and the susceptibility amplitude 
ratio, which have also been the subject of numerous studies over the last several decades. 

2. The determination of pc for site percolation on the square lattice 

Many methods have been developed to estimate percolation thresholds numerically. Some of the more popular 
ones include: 

• Series analysis methods, which involve finding exact statistics for smaller clusters, and were especially impor- 
tant in the earlier years of study ll30l |3TI . 

• The transfer-matrix method |fT2l[r7l . which involves finding exact solutions in finite-width strips, and has been 
recently revived to find very accurate values of pc [29:1 • 

• Finding the value of p where the crossing probability (open system) or wrapping probability (on a torus) equals a 
known amount or becomes independent of the size of the system f2T,^23]. When using the wrapping probability 
on a torus, the convergence of the estimate pdL) is particularly fast: for a square system of size L, one has 

Pc(L) - p, ~ L-'^/* m. 

• Using various real-space renormalization-group theory ideas ll33l[34l . 
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Table 1: Determinations of for site percolation on the square lattice. Numbers in parentheses represent errors in last digit(s). 



year 


author 


method 


Pc 


1960 


Elliott, Heap, Morgan & Rushbrooke yj 


series 


0.48 


1961 


Domb & Sykes [/j 


series 


0.55 


1 C\£. 1 

1961 


Frisch, Sonnenblick, Vyssotsky, Hammersley f3 1 


MC 


0.581(15) 


1963 


Dean ^ 


MC 


0.580(18) 


1964 


Sykes & bssam ||5J 


series 


0.59(1) 


1967 


Dean & Bird 16j| 


MC 


0.591(5) 


1972 


Neal [7] 


MC 


0.593(5) 


1976 


Sykes, Gaunt & Glen ||8J 


series 


0.593(2) 


1976 


Stauiier ||9J 


series 


0.591(1) 


1976 


Leath II 101 


MC 


A COT / 1 /I \ 

0.587(14) 


1978 


T T 1 TT" 1 f) A (T 1 II -1 1 1 

Hoshen, Kopelman & Monberg 111 11 


MC 


0.5927(3) 


1980 


Reynolds, Stanley & Klein lilll 


MC 


0.5931(6) 


1982 


Demda & de Seze 1 12| 


TM 


0.5927(2) 


1982 


• 1' rij. 1 OAiT 1" n 1 oi 

Djordjevic, Stanley & Margolina II 131 


series 


0.5923(7) 


1984 


Gebele L14I1 


MC 


0.59277(5) 


1985 


Rapaport lfT5l 


MC 


0.5927(1) 


1985 


Rosso, Gouyet & Sapoval II16II 


MC 


0.59280(1) 


1985 


Demda & Staufier 11711 


TM 


0.59274(10) 




Z/irr [ ig 1 






1986 


Kertesz lfT9ll 


TM 


0.59273(6) 


1986 


Ziff& Sapoval EO) 


MC 


0.592745(2) 


1988 


Ziff & Stell ED 


MC 


0.5927460(5) 


1989 


Yonezawa, Sakamoto & Hori Il22l 


MC 


0.5930(1) 


1992 


Ziff El 


MC 


0.5927460(5) 


2000 


Newman & Ziff El 


MC 


0.59274621(13) 


2003 


de Oliveira, Nobrega & Stauffer E21 


MC 


0.59274621(33) 


2005 


Deng & Blote ESI 


MC 


0.5927465(4) 


2007 


M. J. Lee E3 


MC 


0.59274603(9) 


2008 


M. J. Lee f28l 


MC 


0.59274598(4) 


2008 


Feng, Deng & Blote ||29ll 


TM/MC 


0.59274605(3) 



• Iterative searching for the value of p where crossing first occurs. The average of that value is used as an estimate 
for/7, Il35l . 

• Adding bonds to a system one at a time until crossing or wrapping first occurs. Here one must convolve with 
the binomial distribution to find the usual fixed-p ("grand canonical") ensemble Il24l |36l [371 . 

• Frontiers or hull- walks in a gradient |fT^|20ll38l[39l , in which the estimate comes from simulations carried out 
in single runs for a given gradient. Vahd in 2d only. 

• Statistics of individual cluster growth with a maximum size cutoff, which eliminates boundary effects. Valid in 
all dimensions, but requires separate runs at different values of p ll40ll4Tl . 

In Table [T] we list the values that have been found for p^. On the average, the precision has improved about one 
digit per decade. Each new digit generally requires about 1000 times more work, both for statistical purposes (a factor 
of 100) and for quantifying finite-size corrections (about another factor of 10); this rate roughly corresponds to the rate 
of increase of computer power (speed and memory) over the years. Two very precise recent works Il28ll29l confirm 
that the value p^ x 0.592746, which was proposed more than 20 years ago 11211 l23l and became a standard after its 
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Table 2: Determinations of SI, oj = DCl = (91/48)f2, and Ai = Sl/cr = (91/36)n. Numbers in parentheses represent errors in last digit(s). 



year 


author 


method 


il 


0) 


A 

Al 


1976 


Gaunt & Sykes |130J 


series 


0.75(5) 


1.42 


1.90 


1978 


Houghton, Reeve & Wallace ||421 


field theory 


0.54-0.68 


0.989-1.28 


1.32-1.71 


1979 


Hoshen et ai. 1143 II 


MC 


0.67(10) 


1 T7 

1.27 


1.69 


1980 


Peai-son |44|| 


conjecture 


64/91-0.703 


1.333 


1.778 


1980 


Nakanishi & Stanley 1145 II 


MC 


0.6 < Q < 1 






1982 


Nienhuis |46l 


field theory 


96/9 Ix! 1.055 


2 


2.667 


1982 


Adler, Moshe & Frivman 113 111 


series 


0.5 


0.95 


1.26 


1983 


A 11 A (T 1 O T~\ a A^B 

Adler, Moshe & Privman 11471 


series 


0.66(7) 


1.25 


1.67 






series for p < pc 


0.49 


0.93 


1.24 




Anarony <x risner ii4o[|4yi 


analytic correction 


cc /Q1 f:r\A 


cc //I 1 1 C 


Jj/jO~i.JJ 


1984 


Margolina et al. ll50ll5T1l 


MC 


0.64(8) 


1.21 


1.62 






series 


0.8(1) 


1.52 


2.02 


1985 


Adler ||49l 


series 


0.63 


1.19 


1.59 


1986 


Rapaport 1521 


MC 


0.71-0.74 






1998 


MacLeod & Jan JSH 


MC 


0.65(5) 


1.23 


1.64 


1999 


Ziff & Babalievski ISll 


MC 


0.77(2) 


1.46 


1.95 


2001 


Tiggemann 1551 


MC 


0.70(2) 


1.33 


1.77 


2003 


Aharony & Asikainen 1561 l57l 


theory, based on l58l 


72/91 


3/2 


2 


2007 


Tiggemann |59| 


MC 


0.73(2) 


1.38 


1.85 


2008 


Kammerer, Hofling, Franosch ll60ll 


MC 


0.77(4) 


1.46 


1.95 


2011 


Ziff 161 1 


theory, based on l62l 


72/91-0.791 


3/2 


2 



inclusion in ||351 . is accurate to all significant figures — and in fact, the next digit is most likely a zero. The average of 
the results of l28]|29l give p, = 0.59274602(4). 

The square lattice site threshold is just one example; thresholds have been studied for scores of systems in various 
dimensions, and many of these results are summarized in the web-page: 

http : / /en . wikipedia. org/wiki/Percolation_threshold. 
Examples of other systems that have been studied extensively include bond percolation on the kagome lattice ||29ll22l 
|63llMllSl|66l and both site and bond percolation on 3d cubic lattices iniE4lHS]|4T]|67l|6i|69ll20ll2ll- 

3. Corrections to scaling 

Many of the methods mentioned above depend upon knowing the behavior of the corrections to the size distribution 
(number of clusters of size s) tisip) at pc, which are expected to be of the form: 

nApc)-As-'(l+Bs-''...) (1) 

where r = 91/48 is universal and is known exactly in 2d. Q is also expected to be universal and has been studied 
in numerous works, usually in the context of determining another quantity such as p^. The resulting values are 
summarized in Table |2] The methods used include series analysis and Monte Carlo simulation on various systems, as 
well as predictions from theory. There have been two recent conjectures that the value of Q is exactly 72/91, the first 
ll56l based upon den Nijs' early result on the corrections to the correlation function [381, and the second |61 1 based 
upon Cardy's more recent result for the crossing probability on an annulus |62|. Recent numerical results |61 1 support 
this value, as reproduced here in Fig.[T] As seen in Table|2] previous estimates and predictions have ranged from 0.49 
to 1.055, but more recent measurements have been closer to the predicted 72/91 ~ 0.791. The exponents Q., cj = DO., 
and Al — Q./cr (where D = 91/48 and cr = 36/91) figure in many problems in percolation. 
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-0.85 

-0.9 
-0.95 



: 


site-SQ 

site-TR 
— »— bond-SQ 

bond-TR 

-72/91 

theory (12) 


A 





2 4 6 8 10 

ln(s/So) 



Figure 1: Plot of n(est) = - log2[(Cj-Cj/2)/(Cs/2 -Cs/4)] where C, = s'^ ^P>jandP>j = 2,!'>,! is the probability of growing a cluster of size 
> s. Data from refs (53|6T]. "Theory" refers to Eq. (12) of gl]. The non-universal metric factor sq equals 0.25 (site-square), 0.13 (site-triangular), 
0.25 (bond-square) and 0.5 (bond-tiiangular). For large s, we see good evidence for the prediction n = 72/91 I56II61I . 



4. Amplitude ratio 

Universality of the scaling function in percolation leads to the prediction that the ratio of the second-moment (or 
the "susceptibility") for equal intervals below and above the threshold is a universal constant, written as C^/C" (where 
"below" is + because of the mapping to the Potts model and corresponds to being above the transition temperature), 
also written as F/F' and F /F+. Finding this ratio has been a notoriously difficult problem in percolation ll93l . as can 
be seen in the summary in Table [3] Recently, three fundamental approaches — series analysis |^T|, Monte Carlo ||9T1 . 
and theory |92| have all converged to a remarkably consistent value of about 161.5(2.0). Note that in some earUer 
work it was speculated that models of continuum percolation may have a different value for this ratio Il75ll80ll83l , 
but more recent work gave results consistent with this value f88l, confirming the expected universality. Likewise, it 
appears that some kinetic percolation-like models, such as kinetic gelation, also have a similar susceptibility amplitude 
ratio. 



5. Conclusions 

We have seen that it has taken decades of work to find consistent values of Q and /C for 2d percolation, 
and to find a very precise value of pc for site percolation on the square lattice. Many other thresholds in various 
dimensionalities are known to relatively low accuracy, and Q. and /C^ are not known to high accuracy in three 
and higher dimensions, leaving much work for future studies. Other amplitude ratios, such as those relating to the 
correlation length, have not been studied to high accuracy, and controversies remain |93 1. Another area for future work 
is precise measurements and characterization of the scaling function in all dimensionalities, not just of the amplitude 
ratios that follow from it. 
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Table 3 : Determinations of tlie amplitude ratio C'*' jC for 2d percolation. 



year 


author 


system, method 




1976 


Sykes, Gaunt, Glen [SJ 


lattice, series (12-20 order) 


1.3-2.0 


1976 


Stauffer ||9l 


lattice, series analysis 


^ 100 


1978 


XT 1 'I'D PI i 1 Il'T'^l 

Nakanishi & Stanley II72I 


lattice, MC 


25(10) 


1978 


Wolff & Stauffer |I221 


lattice, series, fit fiz) to Gaussian 


180(36) 


1979 


Hoshen, Staufier, Bishop, Harrison, Quinn II43I 


lattice, MC 


196(40) 


1980 


Nakanishi & Stanley f431 


lattice, MC (reanalyze) 


219(25) 


1980 


Ahai-ony |74| 


e - b - d expansion 


3.617 


1981 


GawUnsky & Stanley ||75] 


overlapping disks, MC 


50(26) 


1985 


Rushton, Family & Herrmann II76I 


additive polymerization, MC 


140(45) 


1987 


Meir |77| 


lattice, series 


210(10) 


1987 


Kim, Herrmann, Landau llTSjl 


continuum model, MC 


14(10) 


1987 


Nakanishi |179I 


AB percolation, MC 


139(24) 


1988 


Balberg ||80l 


widthless sticks, MC 


=s 3 


1988 


Ottavi |8r| 


approx. theory (Gaussian fit) 


193.9 


19a9 


Corsten, Jan & Jerrard llozl 


lattice, MC 


/j(-l-40, -Ij) 


1990 


S. B. Lee & Torquato ES 


penetrable cone, shell 


1050(32) 


1990 


S. B. Lee il 


disks, MC 


192(20) 


1991 


Hund f85ll 


random contour model, MC 


* 200 


1993 


Zhang & De'Bell Hg) 


Penrose quasi-lattice, series 


310(60) 


1995 


Conway & Guttmann IISTll 


lattice, series (26-33 order) 


45(-h20,-10) 


1996 


S. B. Lee ^ 


penetrable cone, shell, disks 


175(50) 


1997 


S. B. Lee & Jeon in 


kinetic gelation, MC 


170(20) 


1998 


Delfino & Cardy H 


theory, extrapolate Potts to ^ = 1 


74.2 


2006 


Jensen & Ziff ^ 


lattice, MC 


163(2) 


2006 


Jensen & Ziff ED 


lattice, series 


162(3) 


2010 


Delfino, Viti & Cardy USD 


theory, evaluate Potts at = 1 


160.2 
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